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1.0  SUMMARY 


A Terr  adynwaic  Inpact  and  Penetration  Parametric  study  has  been  conduc- 
ted in  which  the  three  principal  parameters,  i.e.,  impact  conditions,  earth 
penetrator  configuration,  *nd  target  media  characteristics,  were  varied  to 
establish  their  effect  on  the  penetration  events  in  question.  The  results  of 
the  study  are  in  the  form  of  trajectory  and  loading  environment  histories 
plus  graphs  which  summarize  performance  trends.  Performance  trends  refer  to 
selected  event  characteristics  such  as: 

• rigid  body  accelerations, 

• penetration  depths, 

• local  pressures  and  shears,  and 

• ricochet  characteristics. 

A brief  description  of  the  computer  codes  used  to  conduct  this  study  is 
provided  in  Section  3.0. 

A detailed  discussion  of  the  study  and  its  results  is  provided  in 
Sections  4 and  5.  In  Section  5.0,  summary  curves  and  tables  are  provided 
which  convey  the  basic  trends  and  results  of  the  study.  The  detailed  tra- 
jectory histories  will  be  provided  under  separate  cover  to  DNA. 


2.0  INTRODUCTION 


The  main  objectives  of  this  study  are: 

1.  Determine  if  the  differential  force  lnw  was  capable  of  correctly 
accounting  for  variations  in  target  media. 

2.  Establish,  parametrically,  a conceptual  understanding  of  the  loading 
environments  imposed  on  an  earth  penetrator  from  impact  through 
detonation. 

The  basic  mission  of  an  EP  includes: 

• delivery  to  the  surface  of  the  target  with  a given  set  of  impact 
conditions,  i.e.,  velocity,  obliquity,  and  angle  of  attack. 

• the  impact  phase  during  which  peak  structural  design  loads  will 
usually  be  experienced. 

• the  penetration  phase. 

• detonation. 

In  order  to  evaluate  and  develop  an  EP  weapon  system,  the  TIfitP  per- 
formance characteristics  must  be  well  defined.  TI&P  characteristics  arc 
required  for  two  important  reasons.  These  are: 

1.  Impact  Loading  Environments 

These  are  necessary  to  EP  structural  design,  component  survival,  aud 
functioning. 

2.  EP  Terradynamic  Performance 

This  information  is  required  for  overall  system  performance  and 
effectiveness  evaluation. 

With  a complete  parametric  data  base  (of  which  this  study  is  only  a 
part)  consisting  of  the  two  performance  characteristics,  other  aspects  associ- 
ated with  the  development  of  an  EP  weapon  system  can  be  addresses.  These 
other  aspects  include: 

• payload  and  component  designation.  • target  selection. 

• EP  packaging  studies.  • target  vulnerability. 

• EP  structural  design.  • EP  effectiveness. 

• component  design. 
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• collateral  damage. 

• delivery  system  considerations. 

For  this  study,  the  TI&P  characteristics  were  established  through  the 
2D  and  3D  Terradynamic  Impact  and  Penetration  Simulator  Computer  Codes  de- 
veloped at  Avco  Systems  Division.  These  codes  are  described  in  detail  in 
Section  3.0. 

The  parametric  study  was  divided  into  six  parts,  each  part  being  selec- 
ted for  the  purpose  of  investigating  certain  EP  performance  phenomena.  The 
'r  arts  include  the  following  series. 

Series  1 

To  determine  the  sensitivity  of  media  variations  which  are  represented 
by  resistance- to-penetration  phenomena  used  in  the  2D  and  3D  codes. 

Series  2 

To  determine  the  effects  of  ballistic  parameter  (i.e. , projectile  weight/ 
cross  section  area)  and  velocity  on  TI&P  performance, 

Series  3 

To  determine  the  effects  of  obliquity  and  angle  of  attack,  on  TI&P 
performance. 

Series  4 

To  determine  the  effects  of  L/D  (i.e. , EP  length/diameter)  on  terra- 
dynamic  stability. 

Series  5 

To  investigate  the  relationships  between  basic  media  engineering  proper- 
ties and  resistance-to-penetration  phenomena,  and  to  predict  the  per- 
formance of  a specific  impact  event  to  be  conducted  at  the  Tonapa  test: 
site  in  New  Mexico. 

Series  6 

To  determine  the  ricochet  performance  characteristics  of  the  Mk  82  and 
Mk  84  bomb  warheads. 

The  specific  matrix  of  runs  for  these  series  are  provided  in  Section  4.0 
along  with  example  results  from  each  series.  The  complete  parametric  study 
results  are  provided  under  separate  cover  to  DNA.  The  parametric  study  re- 
sults are  summarized  and  discussed  in  Section  5.0. 
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3 . 0 IMPACT  AND  PENETRAT ION  SIMULATOR 


Avco's  approach  to  the  analytical  simulation  (or  prediction)  of  a terra- 
dynamic  impact  and  penetration  event  involves  the  use  of  a differential  force 
law.  The  derived  forces  or  pressures  (both  normal  and  tangential)  are  inte- 
grated over  the  surface  of  any  specified  body  of  revolution  to  establish  re- 
sulting loads.  The  equations  of  motion  are  written  (1,  2,  or  3 dimensional'! 
and  are  solved  to  obtain  rigid  bcdy  acceleration,  velocity,  and  trajectory 
histories.  The  approach  is  basically  the  same  as  used  in  present  aerodynamic 
trajectory  analysis  problems,  except  for  the  specific  form  of  the  force  law, 
and  the  discontinuous  nature  of  a typical  impacted  media. 

Avco  has  developed  these  terrao’ynamic  trajectory  programs  for  considera- 
tion of  the  unique  characteristics  of  typical  target  media  and  to  incorporate 
what  we  call  Avco's  General  Resistance-to-Penetration  Force  Law.  This  force 
law  simply  states  the  assumed  form  of  several  possible  resistance-to- 
penetration  phenomena  that  can  occur  during  a terradynamic  impact  event. 

These  phenomena  include: 

• basic  structural  resistance  to  penetration. 

• equivalent  fluid  flow. 

• compressibility  effects. 

• surface  effects. 

The  analytical  form  of  this  force  law  is: 


Normal  Pressure 


Shear  Pressure 
where: 


ur  0 0 

— - = Tj  + 1/2  p CN  sin^  £+  k Ve— a r'  (t-  r) 
dA  1 ’ 

dF 

= fcJj  + 1/2  pV2  Cr  sin  £ c'-'s  £ + fck  V(~°  1 (t  — r) 

dA 


i/  * Basic  structural  resistance  to  penetration 
p «=  Media  density 

V « Local  absolute  velocity  vector 

CN  ■=  Normal  equivalent  flow  coefficient 

C = Local  Incidence  angle 


'Sp' 


k **  Relative  acoustic  impedance  between  impacted  media  and  penetrator 

a - Exponential  decay  factor 

(t-r)  - Mathematical  control  on  the  timing  of  compressibility  effects 

fc  = Coefficient  of  friction 

Cr  = Shear  equivalent  flow  coefficient. 

The  level  of  influence  of  these  phenomena  for  any  specific  impact  event  is 
controlled  by  the  values  of  the  functional  arguments  which  can  vary  with 
respect  to  time,  velocity,  or  displacement. 

Over  the  past  several  years,  Avco  has  obtained  values  of  these  co- 
efficients, empirically,  for  many  representative  media  through  laboratory  and 
ballistic  tests.  The  procedures  used  are  equivalent  to  obtaining  the  lift 
drag  coefficients  of  models  in  wind  tunnels  for  the  aerodynamic  problem.  The 
pertinent  laboratory  tests  included: 

• moisture  content  analysis. 

• sieve  analysis. 

• wave  velocity  measurements  in  the  media. 

• static  resistance  to  penetration  tests. 

• penetrometer  tests. 

The  ballistic  tests  included: 

• firing  models  and  scale  projectiles  into  instrumented  media  samples 
at  varying  impact  velocities  and  impact  angles,  and  measuring 
velocity  and  displacement  (angular  and  translational)  histories, 

• firing  instrumented  full  scale  projectiles  into  hard  (concrete) 
targets  and  measuring  acceleration  and  structural  strain  histories, 
directly.  These  data  were  obtained  through  telemetry  and  recorded 
on  high  speed  magnetic  tape  systems. 

• performing  reverse  ballistic  testing  in  which  the  target  media  was 
fired  at  instrumented  EP's  and  distributed  strain  histories  were 
obtained  by  conventional  "hard  lined"  instrumentation  systems. 

The  laboratory  tests  specifically  yield  the  structural  resistance-to- 
penetration  terms,  the  degree  of  influence  of  compressibility  effects,  and  a 
first  estimate  of  the  influence  of  the  equivalent  flow  terms. 


Through  data  reduction  and  analysis  of  the  ballistic  tests,  specific 
dynamic  values  of  r]  , fc»  CN»  and  Cr  have  been  obtained.  This  procedure  has 
allowed  Avco  to  characterize  many  target  media  and  to  simulate  the  behavior 
of  many  full  scale  projectiles  impacting  these  media.  Accurate  simulations 
of  the  terradynamic  performance  of  several  rocket  configurations,  artillery 
rounds,  EP's,  etc.,  relative  to  trajectory,  loads  and  critical  ricochet  angle 
have  been  made  and  agree  well  with  actual  field  tests  over  the  past  several 
years. 

This  approach  is  relatively  straightforward  and,  similar  to  aerodynamic 
loads  and  trajectory  analysis  techniques,  gives  reasonably  accurate  engineer- 
ing solutions  to  complex  problems  which  otherwise  would  remain  unsolved. 

These  codes  were  used  exclusively  to  conduct  the  parametric  study  described 
in  this  report.  This  technique  which  represents  the  only  capability  for 
efficiently  analyzing  3-D  penetration  problems  is  not  without  its  limitations. 
Compared  with  the  finite-difference  technique  now  available,  the  force  law 
approach  cannot  produce  as  much  detail  in  a given  calculation.  But,  it  is 
important  to  note  that  the  current  finite-difference  calculations  are  limited 
to  axisymetric  calculations  and  the  cost  is  one  to  two  orders  of  magnitude 
greater  than  a corresponding  force  law  calculation. 

The  limitations  of  the  present  force  law  technique  are  associated  with 
the  simplification  regarding  media  response,  distributed  pressures  and  shears 
are  based  upon  the  character  of  the  formulations  presented  above,  the  direction 
and  level  of  the  assumed  total  velocity  field  (which  is  based  solely  upon 
rigid  body  motion  of  the  projectile)  and  whether  a given  location  on  the 
projectile  is  in  the  shadow  of  the  assumed  parallel  flow  field  or  not.  The 
finite-difference  approach  models  the  response  of  the  media  and  therefore 
computes  such  phenomena  as  separation,  compaction,  etc. 
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4.0  EARTH  PENETRATDR  PARAMETRIC  STUDY 


As  discussed  in  Section  2.0,  the  terradynamic  parametric  study  was 
divided  into  six  series,  the  purpose  of  and  the  specific  conditions  of  which 
are  listed  in  Tables  1-6  and  7-12. 

Two  sets  of  tables  are  provided:  the  first  set  being  the  initial 
proposed  parametric  study  and  the  second  (identified  by  "M"  for  modified) 
set  consisting  of  the  deletions  and  additions  to  the  final  matrix  of  runs 
to  be  conducted.  Wherever  possible,  the  original  numbering  system  was 
retained  to  identify  specific  cases. 

The  principal  target  media  for  which  this  study  was  conducted  was  a 
particularly  hard  (i.e.,  high  resistance- to-penetration)  concrete.  Concrete 
resistance-to-penetration  can  vary  by  as  much  as  200  percent  depending  on 
its  age,  basic  unconfined  compressive  strength,  and  the  type  and  size  of  the 
rock  aggregate  employed  in  the  mix.  The  resistance-to-penetration  parameter 
used  in  the  general  resistance-to-penetration  force  law  (described  in  Section 
3)  is  'V  and  thi3  parameter  was  assigned  a value  of  50,000  psi  for  the 
majority  of  the  study.  To  relate  this  value  to  some  equivalent  homogeneous 
media  of ; known  physical  characteristics,  some  conversion  technique  is  re- 
quired. The  "if"  term  is  representative  of  the  local  applied  pressure  exerted 
on  the  EP  which  causes  media  failure  neglecting  inertial  effects.  Any  proven 
model  which  has  been  formulated  to  calculate  this  pressure  can  be  used  to 
relate  the  value  of  the  "if"  term  to  media  physical  properties.  Such  a method 
is  the  so-called  "Cavity  Expansion  Theory"  (first  used  for  the  penetration 
problem  by  J.  N.  Goodier  of  Stanford  Research  Institute).  The  basic  terms 
in  the  "Cavity  Expansion  Theory"  are  similar  to  Avco's  force  law  and  include: 

• basic  media  resistance-to-penetration. 

• inertial  effects. 

• apparent  mass  effects. 

There  is,  therefore,  a one  to  one  correspondence  between  the  two  models  in  at 
least  the  first  terms. 
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TABLE  2 


SERIES  NO.  2 - PEEL. 
VARIATIONS  IN  W/A  AND  VELOCITY 


Problem 

Number 

Type 

2D  or  3D 

Projectile  Data 

Velocity 

(ft/sec) 

Obliquity 

y 

(deg) 

Attack 

Angle 

u 

(deg) 

Material  Property 
Standard  Concrete 
5000  psi 

W/A 

(psi) 

D 

(in) 

Nose 

Remarks 

CRH 

Weight 

(lb) 

2-1 

2D 

6.5 

9.25 

500 

0 

0 

Std 

2-2 

2D 

6.5 

9.25 

1000 

0 

0 

Std 

2-3 

3D 

1S|H 

6.5 

9.25 

1500 

0 

0 

Std 

2-4 

2D 

BfW 

6.5 

9.25 

2000 

0 

0 

Std 

2-5 

2D 

12.1 

6.5 

9.25 

| 

2500 

0 

0 

Std 

2-6 

3D 

12.1 

6.5 

9.25 

400 

3000 

0 

0 

Std 

2-7 

2D 

6.0 

6.5 

9.25 

200 

1000 

0 

0 

Std 

2-8 

2D 

Ha 

6.5 

9.25 

200 

2000 

0 

0 

Std 

2-9 

2D 

m 

6.5 

9.25 

200 

3000 

0 

0 

Std 

2-10 

2D 

18.0 

6.5 

9.25 

600 

1000 

0 

0 

Std 

2-11 

2D 

18.0 

6.5 

9.25 

600 

2000 

0 

0 

Std 

2-12 

2D 

18.0 

6.5 

9.25 

600 

3000 

0 

0 

Std 

2-13 

2D 

24.0 

6.5 

9.25 

800 

1000 

0 

0 

Std 

2-14 

2D 

24.0 

6.5 

9.25 

800 

2000 

0 

0 

Std 

2-13 

2D 

24.0 

6.5 

9.25 

800 

3000 

0 

0 

Std 

2-16 

2D 

30.0 

6.5 

9.25 

1000 

1000 

0 

0 

Std 

2-17 

2D 

30.0 

6.5 

9.25 

1000 

2000 

0 

0 

Std 

2-18 

2D 

30.0 

6.5 

9.25 

1000 

3000 

0 

0 

Std 

2-19 

2D 

31.8 

4.0 

9.25 

400 

1000 

0 

c 

Std 

2-20 

2D 

31.8 

4.0 

9.25 

400 

2000 

0 

0 

Std 

2-21 

2D 

31.8 

4.0 

9.25 

400 

3000 

0 

0 

Std 

2-22 

2D 

8.0 

8.0 

9.25 

400 

1000 

0 

0 

Std 

2-23 

2D 

8.0 

8.0 

9.25 

400 

2000 

0 

0 

Std 

2-24 

2D 

8.0 

8.0 

9.25 

400 

3000 

0 

0 

Std 

TABLE  3 


t SERIES  NO.  3 - PREL. 

f I VARIATIONS  IN  OBLIQUITY  AND  ANGLE  OF  ATTACK 

'/  ’ ! 


TABLE  3 (Concl'd) 


SERIES  NO.  3 


Tyne 

Problem  2D 
Number  3D 


Projectile  Data 


Noe a Remark* 


CRN 


9.25 

With  - Tip 
With  - Tip 
With  • Tip 
With  • Tip 
With  ■ Tip 
With  ■ Tip 


SLA 

Canadian 

400 

400 

400 

400 

400 

400 

400 

400 

Same 

L/D  Wt 

*8 

su 

Same 

L/D  Ut 

08 

su 

Same 

L/D  Wt 

08 

su 

Same 

L/D  Wt 

08 

su 

Same 

L/D  Wt 

08 

su 

Same 

L/D  Ut 

08 

su 

Somo 

L/D  Wt 

08 

su 

Same 

L/D  Wt 

08 

su 

Same 

L/D  Wt 

08 

su 

Same 

L/D  Ut 

08 

su 

Same 

L/D  Ut 

48 

su 

Some 

L/D  Wt 

08 

su 

Same 

L/D  Wt 

08 

su 

Some 

L/D  Wt 

08 

su 

Same 

L/U  Wt 

08 

su 

SLA 

Canadian  400  lb 

SLA 

Canadian  400  Xb 

SLA 

Canadian  400  lb 

su 

Canadian  400  lb 

SLA 

Canadian  400  lb 

su 

Canadian  400  lb 

aoo 

800 

800 

800 

800 

Velocity 

(ft/aec) 


Obliquity 

u 

(deg) 


Attack 

Angle 

a 

(deg) 

material  Property 
Standard  Concrete 
uc  - 5000  psi 

0 

Std 

0 

Ltd 

0 

Std 

0 

Std 

0 

Std 

0 

Std 

0 

Std 

0 

Std 

0 

Std 

+2 

Std 

-2 

Std 

0 

Std 

+2 

Std 

-2 

Std 

0 

Std 

0 

Std 

0 

Std 

0 

Std 

0 

Std 

+2 

Std 

-2 

Std 

0 

Std 

0 

Std 

C 

Std 

0 

Std 

0 

Std 

0 

Std 

0 

Std 

0 

Std 

P2 

Std 

+2 

Std 

+2 

Std 

+2 

Std 

+2 

Std 

1 

5-14 
VI 5 
5-16 
5-17 
5-18 
5-19 
5-20 
5-21 
5-22 
5-23 
5-24 
5-25 
5-26 
5-27 
5-28 
5-29 
5-30 
5-31 
5-32 
5-33 
5-34 
5-35 
5-36 
5-37 
5-38 
5-39 
5-40 


Problem 

Number 

Type 

2D 

3D 

Projectile  Date 

Velocity 

(ft/eec) 

Obliquity 

(deg) 

Attack 

Angle 

(deg) 

Material  Property 

U/A 

(pel) 

D 

(in) 

Noee 

Remark* 

CRH 

SLA 

Canadian 

5-1 

2D 

12.1 

6.5 

9.25 

400 

1500 

0 

0 

oc  - 5000  pel  concrete 

5-2 

2D 

12.1 

6.5 

9.25 

400 

1500 

20 

+2 

uc  - 5000  pel  concrete 

5-3 

2D 

12.1 

6.5 

9.25 

400 

2500 

0 

0 

</c  ■ 5000  pel  concrete 

5-4 

2D 

12.1 

6.5 

9.25 

400 

1500 

0 

0 

ac  » 3000  pel  concrete 

5-5 

2D 

12.1 

6.5 

9.25 

400 

1500 

20 

+2 

o,.  - 3000  pet  concrete 

5-6 

2D 

12.1 

6.5 

9.25 

400 

2500 

0 

0 

0C  • 3000  pel  concrete 

5-7 

2D 

12.1 

6.5 

9.25 

400 

1500 

0 

0 

"Atphalt" 

5-8 

2D 

12.1 

6.5 

9.25 

400 

1500 

20 

+2 

"Aaphalt" 

5-9 

2D 

12.1 

6.5 

9.25 

400 

2500 

0 

0 

"Asphalt" 

5-10 

2D 

12.1 

6.5 

9.25 

400 

1500 

0 

0 

"Till" 

5-11 

2D 

12.1 

6.5 

9.25 

400 

1500 

20 

+2 

"Till" 

5-12 

2D 

12.1 

6.5 

9.25 

400 

2500 

0 

0 

"Till" 

5-13 

2D 

*>1 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 


> 


RESERVED  FOR  ATTEMPTED  SIMULATION  OF  SELECTED  ROCK  PENETRATION 
TESTS  TO  BE  CONDUCTED  BY  SLA  IN  SPRING  OF  FY75 


> 


J 


RESERVED  FOR  SPECIAL  ORDER  OF  DR.  BOB  ROUANI,  UESSD 


Ptobles 

Number 

Type 

2D 

3D 

6-1 

2D 

6-2 

2D 

6-3 

2D 

6-4 

2D 

6-5 

2D 

6—6 

2D 

6-7 

2D 

6-8 

2D 

6-9 

2D 

6-10 

2D 

6-11 

2D 

6-12 

2D 

6-13 

2D 

6-14 

2D 

6-15 

2D 

6-i6 

2D 

6-17 

2D 

6-18 

2D 

6-19 

2D 

6-20 

2D 

6-21 

2D 

6-22 

2D 

6-23 

2D 

6-24 

2D 

TABLE  6 

SERIES  NO.  6 - PREL. 
RICOCHET  PREDICTIONS 


Projectile  Data 

W/A 

D 

Nose  Remarks 

(psi)  (An) 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

To 

Be 

Selected 

(ft/sec) 


Obliquity 

>' 

(deg) 


Attack 

Angle 

U 

(deg) 


Material  Property 
Standard  Concrete 
«c  * 5000  psi 


Hwittr ' 


aatecagfcaA iiiikw? ft  eu /inn 


TABLE  S 


SERIES  NO.  2 - FINAL 
VARIATIONS  IN  W/A  AND  VELOCITY 


Problem 

Number 

Typ« 

20  or  30 

2-1 

2D 

2-2 

20 

2-3 

30 

2-4 

2D 

2-5 

20 

2-6 

30 

2-7 

2D 

2-8 

20 

2-9 

2T 

2-10 

20 

2-11 

2-12 

2D 

2D 

2-13 

2D 

2-14 

20 

2-15 

20 

2-16 

20 

2-17 

20 

2-18 

20 

2-19 

2D 

2-20 

20 

2-21 

20 

2-22 

2b 

2-23 

20 

2-24 

2D 

2-25 

20 

2-26 

2D 

2-27 

2D 

2-28 

2D 

2-29 

2i J 

2-30 

20 

2-31 

20 

2-32 

2D 

2-33 

2D 

2-34 

20 

2-35 

2D 

2-36 

2D 

2-37 

2D 

2-38 

2D 

2-39 

2D 

2-40 

2D 

2-41 

20 

2.-42 

2D 

2-43 

21 

Projectile  Data 


Attack 

Obliquity  Angle  Material  Property 


Standard  Concrete 
5000  pai 


Bun 

Nuabor 

2-3-1 

to  2 

-3-10 

2-4-1 

to  2 

-4-5 

2-5-1 

to  2 

-5-5 

2-6-1 

to  M 

-6-25 

2-7-1 

tO  4 

-7-5 

2-6-1 

to  4 

-8-5 

2-9-1 

to 

-9-5 

2-10-1 

to 

2-10-5 

2-11-1 

to 

2-11-5 

2-12-1 

to 

2-12-5 

2-13-1 

to 

2-13-5 

2-14-1 

to 

2-14-5 

2-15-1 

to 

2-15-5 

2-16-1 

to 

2-16-5 

2-17-1 

to 

2-12-5 

2-18-1 

to 

2-18-5 

2-19-1 

to 

2-19-5 

2-20-1 

to 

2-20-5 

2-21-1 

to 

2-21-5 

2-22-1 

to 

2-22-5 

2-23-1 

to 

2-23-5 

2-24-1 

to 

2-24-5 

2-25-1 

to 

2-25-5 

2-26-1 

to 

2-26-5 

N 

I 

►J 

1 

to 

2-27-5 

2-28-1 

to 

2-26-5 

2-29-1 

to 

2-29-5 

2-30-1 

to 

2-30-5 

2-31-1 

to 

2-31-5 

2-32-1 

to 

2-32*5 

2-33-1 

to 

2-33-5 

2-34-1 

to 

2-34-5 

2-35-1 

to 

2-35-5 

2-36-1 

to 

2-36-5 

2-37-1 

to 

2-37-5 

2-38-1 

to 

2-38-5 

2-39-1 

to 

2-39-5 

2-40-1 

to 

2-40-5 

2-41-1 

to 

2-41-5 

2-42-1 

to 

2-42-5 

2-43-1 

to 

2-4  3-5 

Suanary 
Curve  Mo. 


mmmm 


Problem 

Number 


Projectile  Data 


Nose  Remarks 


TABLE  8 (ConclM) 


Attack 

Obliquity  Angle 
Velocity  y a 

(ft/sec)  (deg)  (deg) 


Summary 
Curve  No. 


2-44-1 

to 

2-44-5 

2-45-1 

to 

2-45-5 

2-46-1 

to 

2-46-5 

2-47-1 

to 

2-47-5 

2-48-1 

to 

2-48-5 

2-49-1 

to 

2-49-5 

2-50-1 

to 

2-50-5 

2-51-1 

to 

2-51-5 

2-52-1 

to 

2-52-5 

2-53-1 

to 

2-53-5 

2-54-1 

to 

2-54-5 

2-55-1 

to 

2-55-5 

2-56-1 

to 

2-56-5 

2-57-1 

to 

2-57-5 

2-58-1 

to 

2-58-5 

2-59-1 

to 

2-59-5 

2-60-1 

to 

2-60-5 

H 

1 

*~4 

vO 

> 

C4 

to 

2-61-5 

2-62-1 

to 

2-62-5 

2-63-1 

to 

2-63-5 

2-64-1 

to 

2-64-5 

2-1A-1 

2-2A-1 

to 

2-2A-5 

2-3A-1 

to 

2-3A-5 

2-4A.-1 

to 

2-4A-5 

2**  5A- 1 

to 

2-5A-5 

2-6A-1 

to 

2-6A-5 

2-1B-1 

to 

2- IB- 5 

2-2B-1 

to 

2-2B-5 

2-1B-1 

to 

2-3B-5 

2-4B-1 

to 

2-4B-5 

2-5B-1 

to 

2-5B-5 

2-6B-1 

to 

2-6B-5 

2-1C-1 

to 

2-1C-5 

2-2C-1 

to 

2-2C-5 

2-3C-1 

*0 

2-3C-5 

2-4C-1 

to 

2-4C-5 

2-5C-1 

to 

2-5C-5 

2-6C-1 

to 

2-6C-5 

2-18-1 

to 

2-18-18 

SERIES  NO.  4 - FINAL 
EFFECT  OF  L/D  ON  STABILITY 


P rob lam 
Number 

Typ® 

2D 

3D 

Projectile  Data 

Velocity 

(ft/sec) 

Obliquity 

y 

Meg) 

Attack 

Angle 

a 

Meg) 

Materiel 

Property 

Run 

Number 

Susmary 
Curve  No* 

V/A 

(pai) 

D 

(in) 

Koaa 

Remarks 

CRH 

L/D 

Weight 

(lb) 

4-1 

2r 

Variable 

. 

6.5 

9.25 

10 

400 

1500 

10 

+2 

Till 

4-1-1  to  4-1-5 

55 

4-2 

2D 

Variable 

6.5 

9.25 

8 

320 

1500 

10 

+2 

Till 

4-2-1  to  4-2-5 

56 

4-3 

20 

Variable 

6.5 

9.25 

6 

240 

1500 

10 

+2 

Till 

4-3-1  to  4-3-5 

57 

4-4 

2D 

Variable 

6.5 

9.25 

4 

160 

1500 

10 

+2 

Till 

4-4-1  to  4-4-5 

58 

4-5 

2D 

Vavi  ble 

6.5 

9.25 

10 

400 

2500 

10 

+2 

Till 

4-5-1  to  4-5-5 

59 

4-6 

2D 

Variable 

6.5 

9.25 

8 

320 

2500 

10 

+2 

Till 

4-'-l  to  4-6-5 

4-7 

2D 

Variable 

6.5 

9.25 

6 

240 

2500 

10 

+2 

Till 

4-7-1  to  4-7-5 

4-8 

2D 

Variable 

6.5 

9.25 

4 

160 

2500 

10 

+2 

Till 

4-8-1  to  4-8-5 

4-9 

2D 

Variable 

6.5 

Blunt 

3 

320 

1500 

10 

+2 

Till 

4-9-1  to  4-9-5 

4-10 

2D 

Variable 

6.5 

Blunt 

4 

160 

1500 

10 

+2 

Till 

4-10-1  to  4-10-5 

4-11 

2D 

Variable 

6.5 

Blunt 

8 

320 

2500 

10 

+2 

Till 

4-11-1  to  4-11-5 

4-12 

2D 

Variable 

6.5 

Blunt 

4 

160 

2500 

10 

+2 

Till 

4-12-1  to  4-12-5 

Table  - 11 

SERIES  NO.  4 - FINAL 
MEDIA  VARIATIONS 


Projectile  Date 

Type 

2D 

3D 

Nose 

Remarks 

Obliquity 

y 

(^og) 

Attack 

Angle 

Meg) 

Problem 

Number 

W/D 

(pai) 

(in) 

CRH 

Weight 

(lb) 

Velocity 

(ft/8er.) 

Material 

Property 

Run 

Number 

Summary 
Curve  No. 

5-4 

2D 

12.1 

6.5 

9.25 

400 

1500 

•j 

0 

o ■ 3000  psi 
Concrete 

5-4-1  to  5-4-5 

5-5 

2D 

12.1 

6.5 

9.25 

400 

1500 

20 

+2 

oQ  - 3000  psi 
Concrete 

5-5-1  to  5-5-5 

60 

5-6 

2D 

12.1 

6.5 

9.25 

400 

250P 

0 

0 

oc  ■ 3000  psi 
Concrete 

5-6-1  to  5-6-5 

61 

5-10 

2D 

12.1 

6.5 

9.25 

400 

1500 

0 

0 

Till 

5-10-1  to  5-10-5 

60 

5-11 

ID 

12.1 

6.5 

9.25 

400 

1500 

20 

+2 

Till 

5-11-1  to  5-11-5 

61 

5-12 

2D 

12.1 

6.5 

9.25 

400 

2500 

0 

0 

Till 

5-12-1  to  5-12-5 

TONAPA  TEST  PREDICTIONS 


5-1 A 

2D 

12.1 

6.5 

9.25 

517 

1500 

0 

0 

Till 

5-1A-1  to  5-1A-5 

66 

5-2A 

2D 

12.1 

6.5 

9.25 

517 

2000 

0 

0 

Till 

5-2A-1  to  5-2A-5 

67 

5-3A 

2D 

12.1 

6.5 

9.25 

517 

1640 

0 

0 

Till 

5-3A-1  to  5-3A-5 

25 


wan— 


L\ 1 V «■  Ji'Arc' xz  Vi 


TABLE  12 


Problem 

Number 


Projectile  Data 


Type 

2D  W/A  D 

3D  (pei)  (in) 


DELETED 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

DELETED 

3D 

3D 

3D 

3D 

2D 

2D 

2D 

DELETED 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

2D 

DELETED 
2D  I I 


SERIES  NO.  6 - FINAL 
RICOCHET  PREDICTIONS 


1 Nobs  Ramarke  | 

Mk 

Waight 

(lb) 

Obliquity 

Attcok 

Angle 

Valocity 

y 

a 

(ft/aec) 

(dag) 

(dag) 

Material  Proparty 


arget  Layerad  Madia 
argat  Layared  Madia 
arget  Layered  Madia 
argat  Layered  Madia 
argat  Layered  Madia 
argat  Layarad  Madia 
argat  Layarad  Madia 
argat  Layarad  Madia 
argat  Layarad  Madia 
argat  Layarad  Madia 
'argat  Layarad  Madia 
'argat  Layarad  Madia 
'argat  Layarad  Madia 
>rgat  Layared  Media 


6-2A-1  to  6-2A-3 
6-3A-1  to  6-3A-3 
6-4A-1  to  6-5A-3 
6-5A-1  to  6-5A-3 
6-6A-1  to  6-6A-3 
6-7A-1  to  6-7A-3 
6-8A-1  to  6-8A-3 
6-9A-1  to  6-9A-3 
6-10A-1  to  6-10A-3 
6-UA-l  to  6-11A-3 
6-12A-1  to  6-12A-3 
6-13A-1  to  6-13A-3 
6-14A-1  to  6-14A-3 
6-15A-1  to  6-15A-3 


Target  Layarad  Madia  6-4-1  to  6-4-10 
Target  Layarad  Madia  6-5-1  to  6-5-12 
Target  Layarad  Madia  6-12-1  to  6-12-6 
Target  Layerad  Madia  6-13-1  to  6-13-10 
Target  Layered  Madia  6— IB— 1 to  6-IB-3 
Target  Layarad  Madia  6-2B-1  to  6-2B-3 
Target  Layarad  Madia  6-3B-1  to  6-3B-3 


Layared  I 
Layerad  I 
Layerad  1 
Layered  1 
Layarad  1 
Layarad  1 
Layarad  ] 
Layared 
Layered 
Layered 


6-5B-1  to  6-5B-3 
6-7B-1  to  6-7B-3 
6-8B-1  to  6-8B-3 
6-9B-1  to  6-9B-3 
6-10B-1  to  6-10B-3 
6-12B-1  to  6-12B-3 
6-14B-1  to  6-14B-3 
6-1C-1  to  6-1C-3 
6-2C-1  to  6-2C-3 
6-3C-1  to  6-3C-3 


Target  Layered  Madia  6-5C-1  to  6-5C-3 
Target  Layerad  Media  6-7C-1  to  6-7C-3 
Target  Layered  Media  6-8C-1  to  6-8C-3 
Target  Layared  Media  6-9C-1  to  6-9C-3 
Target  Layerad  Media  6-10C-1  to  6-10C-3 
Target  Layered  Media  6-12C-1  to  6-12C-3 
Target  Layered  Media  6-14C-1  to  6-14C-3 
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The  equation  of  motion  used  to  describe  the  terradynamic  trajectory  of 
a penetrator  based  on  the  Cavity  Expansion  technique  is: 


d2z  nD2  { 2 f d/cAN  /dz\2l) 

M * * ~ r + T p?  r r\~2f 


The  parameters  in  Eq.  (1)  are  given  in  terms  of  target  properties  by  the 
following  relations: 


« 1-5 


"3 

T 


d + «p)S1/3  + 


<A3  *■  — E - exp  ( -3  /3  )J  - — Y£n  8 


+ — rr^Et-  — Et?7 
27  ' 9 c 


5.1 exp  ( -3  P ) 

'P 


^ 2E  3 


Pp  = P0  «P<V  <9> 

where : 

O 

p0  - initial  density  of  target  material,  slugs/ft 
Y - yield  strength  of  target  material,  psf 

E - Young's  modulus  of  elasticity,  corresponding  to  locked  elastic 
region,  psf  ■ 3G 

Et  - strain-hardening  modulus,  corresponding  to  locked  plastic  region, 
psf  “ 3Gt 

•Rohenl,  R„  Waterway*  Experiment  Paper, S-73-68,  June  1073. 
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ej  ■ volumetric  strain  related  to  the  elastic  region  of  the  pressure- 
volumetric  strain  curve 

ep  - volumetric  strain  related  to  plastic  region  of  the  pressure- 
volumetric  strain  curve 
n - 3.1416 

fN  ■ sin  9 + 1/3  cos^  6 
9 ■ nose  angle 


The  equation  of  motion  using  Avco's  force  law  approach  for  one  dimen- 
sional motion  is: 


-M  — - fc)4-*  1/2  pV2  (Cj^j  sin2  £+-►  Cr  sin  £ cos  £)  dA 


(10) 


Equating  similar  terms: 

2/3  pp  fN  02  V2  - 1/2  p V2  (CN  sin2  £ + Cf  sin  £ cos  £) 


(ID 


2/3  pp  0! 


modified  nose  shape 
4-  increased  mass 


(12) 


03  « rj  fc)  (13) 

Referring  to  Eq.  (11)  both  sides  represent  equivalent  fluid  flow  terms.  The 
Cavity  Expansion  Theory  offers  no  direct  or  fundamental  method  of  obtaining 
the  correct  values  of  the  required  coefficients,  consequently,  empirical 
evaluation  of  these  are  necessary  unless  a conversion  process  can  be  derived 
from,  for  example,  trial  finite  difference  calculations. 


The  Cavity  Expansion  term  represented  by  Eq.  (12)  is  representative  of 
an  equivalent  mass  of  fluid  entrapped  by  the  nose  of  the  panetrator  and  de- 
celerating with  it. 


The  conversion  here  is  direct;  however,  in  the  case  of  a large  sharp 
nose  penetrator,  the  term  becomes  negligible. 


In  the  Cavity  Expansion  Theory  sliding  friction  between  the  media  aud 
the  projectile  surface  is  not  represented.  Consequently,  Eq.  (13)  becomes: 

(14) 
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03  = 7? 


The  ">7"  term  used  In  Avco's  force  lav  approach  is  a measurable  quantity 
established  by  conducting  static  penetration  tests.  "03"  is  the  analytically 
derived  static  resistance  to  penetration.  Consequently , "ij"  can  be  obtained 
by  evaluating  the  "^3"  expression. 

To  summarize,  only  two  of  the  penetrability  terms  for  Avco's  force  law 
approach  can  be  derived  from  the  Cavity  Expansion  technique.  These  include 
the  apparent  mass  term,  which  already  has  been  included,  and  the  normal 
pressure  term  associated  with  the  media's  static  resistance-to-penetration. 
The  sliding  or  Coulomb  friction  term  used  in  Avco's  force  law  technique  is 
not  currently  present  in  the  Cavity  Expansion  technique,  and  the  equivalent 
fluid  flow  terms  used  in  both  techniques  must  be  arrived  at  empirically. 

This  conversion  or  equivalency  technique  can  therefore  be  used  to  de- 
scribe the  strength  characteristics  of  the  concrete  used  in  the  parametric 
study. 

Performing  the  indicated  operations  suggeeted  by  Eq.  (14)  and  assuming 
that  Young's  modulus  remains  that  of  typical  concrete,  l.e.,  3 x 106  psi, 
the  equivalent  unconflned  compressive  strength  is  approximately  10,000  psi. 

It  is  apparent  that  if  the  above  procedure  can  be  used  to  arrive  at  and 
describe  an  equivalent  strength  media,  it  can  also  be  used  to  find  the 
approximate  "q"  term  given  the  properties  of  the  media.  This  technique  has 
been  used  in  the  past  to  aid  in  predicting  penetration  events  into  "un- 
characterized" media  and  an  example  of  this  procedure  and  outcome  is 
described  in  Section  5.0  in  which  a prediction  was  made  for  the  Tonapa 
Test  Series  conducted  in  the  Spring  of  1975. 

The  following  sections  describe  the  results  of  the  parametric  studies 
for  Series  1.0  through  6.0. 
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5.0  IMPACT  AND  PENETRATION  PARAMETRIC  STUDY  RESULTS 
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In  this  section  each  series  of  runs  Is  described  and  summarized  to  allow 
the  reader  to  quickly  grasp  the  more  significant  aspects  of  the  particular 
area  of  Investigation.  Under  conditions  where  a more  In-depth  review  Is 
desired,  the  reader  Is  referred  to  the  detailed  time  histories  of  the  event. 
When  specific  runs  are  referred  by  run  number  (l.e.,  3.8,  3.60  etc.)  the 
variables  are  defined  in  the  appropriate  table  of  this  report  and  the  complete 
run  data  is  available  under  separate  cover  through  DNA. 

5.1  SERIES  1 — SENSITIVITY  STUDY 

The  matrix  of  runs  conducted  for  Series  1 is  presented  in  Table  7.  The 
EP  used  in  this  study  is  shown  in  Figure  1.  The  EP  weighs  400  pounds,  is 
60  inches  long  with  a diameter  of  6.5  inches.  The  nose  shape  is  described  by 
the  parameter  "CRH"  which  represents  the  ratio  between  the  tangent  ogive 
radius  and  the  EP  diameter. 

The  nominal  media  for  this  series  is  the  "hard  concrete"  described  in 
Section  4.0  and  variations  from  the  nominal  are  given  in  percent.  In  each  of 
the  series,  there  are  generally  two  types  of  runs  made;  i.e.,  2-D  and  3-D.  All 
the  impact  conditions  are  described  in  a plane  and  because  none  of  the  EP's 
is  spinning,  no  out  of  plane  motion  is  incurred.  Consequently,  the  impact 
event  can  be  described  in  the  2-D  plane  of  symmetry  which  they  are.  The  pur- 
pose of  the  3-D  runs  is  to  provide  three  dimensional  distribution  plots  of 
the  local  pressures  and  shears. 

Examples  of  the  2-D  and  3-D  plots  for  Runs  3.60  and  3.8  are  presented  in 
Figures  2 through  6 and  7 through  16,  respectively.  The  output  format  for 
these  plots  is  described  in  detail  in  Tables  13  and  14. 

There  are  generally  five  separate  graphs  for  each  2-D  run  consisting  of: 

• a pictorial  description  of  the  event.  (See  Figure  2.) 

• the  rigid  body  resulting  loading  environments.  (See  Figure  3.) 

• the  basic  trajectory  data.  (See  Figure  4.) 

• the  rigid  body  axial  accelerations.  (See  Figure  5.) 

• the  rigid  body  lateral  accelerations.  (See  Figure  6.) 
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Figure  1.  Target  projectile  (EP| 
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Figure  6.  Lateral  acceleration 
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Figure  9.  (3D)  Rotational  data 
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Figure  15.  (3D)  Stritt  data 


PURPOSE: 


TABLE  13 


To  give  the  output  format  for  2D  plota. 


GENERAL  FORMAT:  At  the  top  of  each  ^raph  directly  under  the  main  title 

specific  data  on  the  penetrator  characteristics  and  impact 
conditions  are  included. 


Symbol  Characteristics 

R Cylindrical  radius 

L Cylindrical  length 

GAM1  Impact  obliquity  angle  plus  angle  attack 

VEL  Impact  velocity  of  the  penetrator,  while  it  is 

zero  when  the  angle  of  attack  is  not  zero. 

XCG  Center  of  gravity  acceleration 

XU  Height  of  the  penetrator 

1.  Graph  1:  Trajectory  Schematic 


Units 

inches 

inches 

degrees 

inches/sec 


ft/sec^ 

pounds 


I 1 


The  penetrator  axis  lies  in  the  XZ  plane.  Velocity  of  the  penetrator  is 
shown  as  an  arrow  from  the  center  of  gravity  and  point  along  the  axis  of 
the  penetrator.  Alpha  is  the  angle  of  attack  which  is  defined  as  follows: 

a a y - b ' 

where  y is  obliquity  angle  measured  from  the  normal  of  the  target  to  the 
axis  of  the  penetrator  as  shown  below: 


TABLE  13  (Coat'd) 


and  6 ' la  def load  a* 


0' 


tan 


-1 


AX 

AZ 


therefore  a ie  positive  when  y ■»  O' 
a is  negative  when  y < 6' 


2.  Graph  2:  The  Loading  History  of  the  Penetration 


Symbol 

FXP 

FZP 

MCC 


Characteristics  Units 

The  lateral  loads  pounds 

The  axial  loads  pounds 

The  applied  moment  (or  torque)  acting  about  the  in-pounds 
center  of  gravity. 


3.  Graph  3:  The  Trajectory  History  of  the  Penetration 


Symbol 


Characteristics 


X 

Z 

VEL 

ALA 

GAM 


The  lateral  displacement 

The  axial  displacement 

The  velocity  at  the  center  of  gravity 

The  angle  of  attack  amplitude 

The  angle  of  obliquity 

(orientation  of  the  penetrator  measured  from 
the  target  normal) 


Units 

inches 

inches 

inches/sec 

radians 

degrees 


4.  Graph  4;  The  Axial  Acceleration  History  of  the  Penetration 


Symbol  Characteristics  Units 

GZ1  ' * ' The  axial  acceleration  at  point  1 g 
GZ2  ' * ' The  axial  acceleration  at  point  2 g 
GZ3  ' 4 ' The  axial  acceleration  at  point  3 g 


and  the  three  points  are  referred  to  the  positions  on  the  penetrator 
selected  by  the  user.  (For  this  study  at  the  eg). 
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TABLE  13  (Cone I’d) 


Graph  5 
iymbol 


: The  Lateral  Acceleration  History  of  the  Penetration 
Characteristics 

The  lateral  acceleration  at  point  1 g 
The  lateral  acceleration  at  point  2 g 
The  lateral  acceleration  at  point  3 g 
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TABLE  14 


PURPOSE: 

GENERAL  FORMAT: 


To  give  the  output  fonut  for  3D  plots: 


At  the  top  of  tech  graph  directly  under  the  aain  title 
specific  data  on  the  penetrator  characteristics  and  impact 
conditions  are  included. 


1. 


For  further  reference  see  Figure  17. 

Graph  1:  Penetrator  Orientations  in  ZY  Plane  and  XY  Plane 


Velocity  of  the  penetrator  is  not  shown  on  this  graph,  but  the  definition 
of  obliquity  angle  is  a little  different  (by  adding  180°)  from  the  on> 
defined  in  2D  cases. 


2.  Graph  2: 
Symbol 


The  Loading  History  of  the  Penetrator 
Characteristics 


Units 


FX 

• . • 

The  lateral  loads  in  X-dlrection 

pounds 

FY 

• * • 

The  lateral  loads  in  Y-direction 

pounds 

FZ 

• A » 

The  axial  loads  in  Z-direction 

pounds 

MX 

• C£J  1 

The  shear-moment  in  X-direction 

in-pounds 

MY 

* X • 

The  shear-moment  in  Y-direction 

in-pounds 

MZ 

’ ♦ • 

The  bent-moment  in  Z-direction 

in-pounds 

3. 

Graph  3: 

Rotational  History  of  the  Penetrator 

Symbol 

Characteristics 

Units 

YAW  1 . * 

Angular  displaceafent  with  respect  to  local  non- 

degrees 

spinning  frame  of  Zj, 
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Symbol 

Characteristics 

Units  f 

degrees 

i 

PSII  ’ * 

Initial  rotational  displacement  in  yaw  of  the 
flight  path  frame  with  respect  to  inertia  frame. 

GAMI  1 y • 

Initial  rotational  displacement  in  pitch  of  the 
flight  path  frame  with  respect  to  inertia  frame. 

degrees 

I 

PHII  * <f>  • 

Initial  rotational  displacement  of  body  fixed 
frame  with  reapect  to  flight  path  frasie. 

degrees  j 

VELOCITY 

Impact  velocity  of  the  penetrator,  while  it  is 
zero  when  the  angle  of  attack  is  not  zero. 

I 

in/ sec 

TIME 

Time  at  which  the  plot  is  taken. 

sec. 

g|  ■ i ji  m fn  1 1 1 1 [i  r • t~ti — — - 


TABLE  14  (Cont’d) 


Symbol 

Characteristics 

Units 

PCH  ' * 

f 

Angular  displacement  with  respect  to  local  non- 
spinning frame  of  Z2« 

degrees 

ROL  1 * 

f 

Angular  displacement  with  respect  to  local  non- 
spinning frame  of  Z3. 

degrees 

YAR  'Cfl* 

Rate  of  change  of  angular  displacement  with 
respect  to  Zj  along  the  direction  of  Z3. 

rad/ sec 

PCR  ' X 

f 

Rate  of  change  of  angular  displacement  with 
respect  to  Z2  along  the  direction  of  Zj^. 

rad/sec 

ROR  ' ♦ 

f 

Rate  of  change  of  angular  displacement  with 
respect  to  Z3  along  the  direction  of  Z£. 

rad/sec 

NOTE:  All  rotations  were  made  using  right-hand  rule. 

4.  Graph  4 

: Translational  History  of  the  Penetrator 

Symbol 

Characteristics 

Units 

X ' • ' 

The  lateral  displacement  in  X-direction 

inches 

Y ' * ' 

The  lateral  displacement  in  Y-direction 

inches 

z ' A • 

The  axial  displacement  in  Z-direction 

inches 

ALA  ’M' 

The  axial  amplitude  of  angle  of  attack 

degrees 

ALL  ' X 

f 

The  lateral  amplitude  of  angle  of  attack 

degrees 

VEL  ’ ♦ 

f 

The  resultant  translational  velocity 

ir./sec 

5 . Graph  5 

to  Graph  10:  Stress  History  of  the  Penetrator 

Symbol 

Characteristics 

Units 

P 

Axial  shear  at  a local  surface 

psi 

Q 

Lateral  shear  at  a local  surface 

psi 

R 

Normal  pressure  at  a local  surface 

psi 

P(0°)  ' 

Axial  shear  at  0° 

psi 

Q(0°) 

t 

• 

’ Lateral  shear  at  0° 

psi 

R(0°)  ) 

Normal  pressure  at  0° 

psi 

P(90°)' 

| 

Axial  shear  at  90° 

psi 

Q(90°) 

>’  * 

’ Lateral  shear  at  90° 

psi 

R(90°) 

) 

Normal  pressure  at  90° 

psi 

The  3-D  graphical  output  format  is  similar  to  the  2-D  except  that  the  rota- 
tional and  translational  trajectory  data  are  separated  (to  avoid  confusion) 
and  several  additional  graphs  (the  number  depending  upon  the  times  during  the 
trajectory  of  interest)  are  generally  provided  to  define  the  pressure  and  shear 
distribution  on  the  body. 

The  3-D  graphical  displays  include: 

• a pictorial  representation  of  the  impact  and  penetration  event.  (See 
Figure  7.) 

• the  total  rigid  body  applied  loads  and  moments.  (See  Figure  8.) 

(The  center  of  gravity  accelerations  can  be  obtained  through  division 
of  these  quantities  by  the  EP  weight  or  rotational  inertias.) 

• the  rotational  motion  and  rate  histories.  (See  Figure  9.) 

• the  translational  trajectory  histories.  (See  Figure  10.) 

• local  pressure  and  shear  distribution  plots.  (See  Figures  11  through 

! 16.) 

i 

! For  specific  definition  of  these  displayed  quantities,  the  reader  is  referred 

to  Tables  13  and  14. 

The  2-D  output  information  is  readily  understood  as  to  the  manner  of  the 
quantities  in  question.  In  order  to  understand  the  output  quantities  of  a 
3-D  event,  it  is  necessary  to  understand  the  reference  frames  in  which  the 
data  is  provided.  The  3-D  data  is  specified  in  three  separate  reference 
frames. 

• Inertial. 

• Flight  path. 

• Local  non-spinning. 

These  reference  frames  and  the  coordinate  transformation  sequence  are  pre- 
sented and  defined  in  Figures  17  and  18.  Referring  to  these  figures: 

• the  translational  state  vectors  are  specified  in  the  inertial  refer- 
ence frame  X , Y , ar.4  Z . 

• the  rotational  state  vectors  are  Eulerian  rotational  displacements 
and  rates. 

I 
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ROTATIONAL  SEQUENCE 

07&  0 


Yi  Y2Y3  ~ Flight  path  frame 
^1X2X3  ~ Body  fixed  frame  (spinning) 
Z'j  Z2Z3  ~ Local  non -spinning  frame 


Figure  IB.  Local  non-spinning  frames 


2> 


Y2 

SIGMAN 

TAUA 

TAUB 
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• the  resulting  forces  and  moments  are  specified  in  the  flight  path 
frame  (Yj , Y2  , and  Y3  ) . 

• the  pressures  and  shears  are  defined  in  the  local  non-spinning  refer- 
ence frame  (Zj,  Z2  , and  Z3).  (See  Figure  18.) 

In  order  to  obtain  the  body  fixed  rotational  vectors  (&>},  oj2  , and  ) 
which  lie  along  the  Xj , X2  , and  X3  axis  respectively  in  Figure  17,  the  follow- 
ing relationships  are  provided: 

^ cosy 

<o2  = (sin  y ijt  — y sin  <f>/cos</>  )/(sin2  (j>  + cos^) 

= (y  + <u2  sin</>)/cos  <£ 

When  studying  the  stress  distribution  data  (i.e.,  Figures  11  through  16),  it 
should  be  noted  that  each  of  these  graphs  provides  surface  traction  informa- 
tion (lb /in2): 

• at  a specific  time  during  the  penetration  event. 

• at  specific  locations  around  the  periphery  of  the  EP. 

• distributed  along  the  length  of  the  vehicle  (the  nosetip  being  the 
zero  station)  and  after  stations  along  the  abscissa. 

• consisting  of  three  components,  SIGMAN  - normal  to  the  surface, 

TAUB  - perpendicular  to  the  body  longitudinal  axis  in  the  plane  of 
the  surface,  and  TAUA  - also  in  the  plane  of  the  surface  directed  aft. 

Some  discussion  relative  to  the  shape  of  these  surface  tractions  distri- 
butions is  in  order.  The  specific  value  of  the  pressure  or  shear  is  of  course 
the  direct  vectoral  sum  of  the  force  law  terms  described  in  Section  3.0. 

In  the  majority  of  runs  conducted  for  this  study,  these  terms  consisted 
of  the  media  resistance  terms  (i.e.,  r\  and  f crj  ) and  the  equivalent  fluid  flow 
terms  (C^q,  and  Crq  , where  q is  thf  dynamic  pressure).  Because  of  the  shape 
of  these  EP's  (mostly  sharp  nosed)  little  compressibility  effects  were  noted 


f r*m 


and  only  for  high  obliquity  cases  (i.e,,  45°  and  above)  were  surface  effects 
introduced.  As  a result,  the  surface  traction  data  appear,  in  general,  as 
steady-state  conditions.  This  requires  some  explanation.  Avco  theorizes 
that  there  are  three  penetration  regions  for  a normal  impact  event  into  a 
semi-infinite  target  half  space. 

1.  Initial  region  — where  the  presence  of  the  free  surface  significantly 
effects  distributed  surface  tractions. 

2.  The  transition  region  — where  the  free  surface  is  becoming  lens 
influential. 

3.  The  steady  state  region  — where  the  presence  of  free  surfaces  no 
longer  effects  the  surface  tractions. 

Little  real  test  data  is  available  which  quantitatively  describes  the  first 
two  phases  or  the  transition  between  Phases  1 and  3.  In  addition,  these 
phases  vary  significantly  depending  upon  target  media,  EP  nose  configuration, 
and  velocity.  Fortunately,  the  resulting  (integrated)  loads  experienced  by 
an  EF  during  this  phase  of  the  vent  do  not  appear  to  alter  subsequent  steady 
state  loading  environments  or  penetration  performance  significantly.  For  the 
majority  of  runs  therefore,  the  steady-state  (i.e.,  third  phase)  force  law 
theory  was  used.  This  accounts  for  the  somewhat  square  appearance  of  the  sur- 
fact  tractions  distributions.  Recently,  finite  difference  calculations  of 
the  initial  phase  of  the  penetration  event  have  been  made  and  these  surface 
tractions  are  more  triangular  in  shape.  In  this  parametric  studv,  where  the 
initial  phase  of  the  penetration  event  significantly  effects  subsequent  pene- 
tration performance,  the  initial  phase  theory  is  used.  The  effect  of  this  is 
to  produce  surface  traction  distributions  more  similar  to  those  predicted  by 
the  finite  difference  techniques. 

The  detailed  trajectory  and  loads  history  data  for  the  Series  1,  Sensi- 
tivity Study,  was  indicated  in  Runs  1.1-1  through  1.20.5  as  defined  in 
Table  7.  These  data  have  been  analyzed,  reduced,  and  summarized  in  Figures  19 
through  24  of  this  report. 


Velocity  = 1500  ft/*ec  Angle  of  attack  = 0 degrees 
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Projectile:  W/A  = 12.1  psi, 
CRH  = 9.25.  W 
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Figures  19  and  20  show  penetration  depth  effects  caused  by  changes  in 
the  Avco's  Force  Law  coefficients  (i.e.,CN  , Cr  , fc  and  jj  ) and  angle  of 
attack  where: 

CN  “ Normal  pressure  coefficient  due  to  flow  effects. 

Cr  - Shear  pressure  coefficient  due  to  flow  effects. 

t]  « Media  structural  resistance  to  penetration. 

fc  ■>  Media  coefficient  of  friction  (used  as  f <-*»/)• 

Only  two  variations  of  these  parameters  were  u ed,  (i.e.,  0.0  and  501 
change  from  the  nominal  value) , consequently  only  a straight  line  can  be 
drawn  connecting  the  points.  It  is  apparent  from  these  figure  j,  that  " » " has 
th*  most  significant  effect  on  penetration  depth  achieved.  The  change  in 
depth  caused  by  the  3 degree  angle  of  attack  was  less  than  10  percent  as  shown 
in  Figure  20. 

The  variation  in  resulting  rigid  body  g's  for  the  normal  impact  situation 
caused  by  these  parameter  changes  is  presented  in  Figures  21  through  24.  The 
peak  axial  nominal  value  expected  for  this  impact  situation  would  be  on  the 
order  of  8000  g's.  The  maximum  variation  consists  of  an  increase  to  11,000 
g's,  again  because  of  the  'Y'  change. 

As  the  angle  of  attack  varies  from  zero  to  three  degrees,  the  axial  loads 
will  Increase  because  of  the  addition  of  friction  and  shear  along  the  cylin- 
drical portion  of  the  EP.  (See  Figure  22.)  Side  loads  and  moments  are  also 
generated  as  indicated  in  Figures  23  and  24.  The  variation  of  axial  loading 
environments  caused  by  changes  in  obliquity  is  included  in  Figures  21  and  22. 
Changes  in  the  lateral  loade  and  moments  caused  by  an  obliquity  of  20  degrees 
is  negligible  compared  to  the  effects  of  angle  of  attack.  The  minor  differences 
can  be  noted  by  referring  to  detailed  trajectory  history  results  presented  in 
the  data  of  Volume  I.  The  most  sensitive  parameter  for  these  cases  is  again 
The  largest  variation  occurs  for  the  axial  loading,  again  caused  by  *Y'. 

Less  change  occurs  for  t’ne  lateral  load  and  almost  no  change  occurs  for 
the  resulting  moment.  The  reason  for  this  is  because  in  the  case  of  the  axial 
loads,  each  component  of  the  increasing  load  (i.e.,  normal  and  frictional 
component)  add  thereby  generating  a relatively  large  overall  increase,  in 
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the  case  of  the  lateral  load  and  the  moment,  these  components  subtract  and 
result  in  a smaller  ini  rease. 

Cases  1-2,  1-17,  and  1-18  were  selected  from  Series  1 for  3-D  runs  to 
determine  the  surface  tractions  data.  Peak  values  of  pressures  and  shears 
were  extracted  from  the  basic  data  and  are  provided  in  Table  15.  The  varia- 
tions in  the  peak  levels  are  not  significant.  What  is  significant  follows 
from  the  integration  of  tnese  distributions  which  causes  significant  lateral 
loads  and  moments  to  be  applied  to  the  EP  structure, 

5.2  SERIES  2 - VARIATIONS  IN  W/A  AND  VELOCITY 

The  purpose  of  parametric  simulations  defined  by  Series  2 is  to  establish 
the  effects  on  performance  and  loading  environments  caused  by  variations  in 
"W/A"  (EP  weight  to  cross  sectional  area  ratio,  a ballistic  parameter),  and 
impact  velocity.  The  runs  selected  to  demonstrate  these  effects  are  presented 
in  Table  8.  Two  other  variables  are  included  in  the  study.  These  are  EP 
diameter  and  the  parameter  "CRH",  "CRH"  is  defined  as  the  ratio  between  the 
ogive  radius  and  the  cylinder  diameter  of  the  EP.  "CRH"  therefore  is  a 
measure  of  nose  sleekness.  Series  2 consisted  of  Runs  2.1.1  through  2.64.5 
as  defined  in  Table  8.  The  data  summary  curves  for  Series  2 are  presented  in 
Figures  25  through  36. 

For  a "CRH"  of  9.25  and  an  E?  diameter  of  6.5  inches  the  relationship 
between  depth  of  penetration  and  "W/A"  for  velocity  of  1000,  2000,  and  3000 
ft/sec  is  shown  in  Figure  25.  The  peak  axial  rigid  body  "g"  loads  for  these 
same  variations  are  presented  in  Figure  26.  The  only  significant  trend  noted 
foi  this  series  of  runs  is  related  to  penetration  performance  trends  as  in- 
fluenced by  W/A  as  opposed  to  velocity.  As  a result  of  the  loading  theory 
(i.e.,  being  relatively  constant  which  in  Avco's  view  is  characteristic  of 
penetration  events  into  hard  targets),  increasing  velocity  is  a much  more 
effective  way  to  improve  penetration  performance  than  equivalent  linear 
increases  in  W/A.  At  the  same  time,  .ie  rigid  body  loading  environments,  over 
the  velocity  range  of  interest  are  not  significantly  altered  by  increases 
in  velocity. 


65 


Ir.g'arl 
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Figure  26.  Axial  acceleration:  CRH  - 0.25 
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Figure  28.  Penetration  performance:  CRH  *=  2.25 


Projectile:  CRH  * 6.Q,  Oiamster  * 6.5"  Target:  STD  5000  pti  concrete 


Figure  29.  Penetration  performance:  CRH  * 6.0 


Projectile:  CRH  = 0.5,  Diameter  «=  6.5"  Target:  STC  5000  psi 
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Figure  35.  Penetration  performance:  blunt  nose 
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Figures  27,  28,  and  29  show  how  performance  changes  as  a function  of 
"CRH"  as  all  other  par.-uneters  are  held  constant.  There  is  a crossover  in  the 
effects  of  "CRH''.  This  is  not  a very  significant  anomaly,  but  the  data  does 
imply  that  as  "CRH"  varies  from  9.25  to  6.0  to  2.5  to  0.5  the  penetration 
performance  at  first  improves  a little  and  then  pr  ceeds  to  get  worse  as  the 
EP  becomes  fairljr  blunt.  This  effect  is  caused  by  the  analytical  formulation 
of  Avco's  force  xaw  and  specifically  by  how  bluntness  (as  controlled  by  "CRH") 
is  treated  in  each  term  of  1 .e  *orce  law;  i.e.: 

The  Impact  Shock  Term  ■ f (sin  0)  , i.e. , is  a function  of  sin  0 

Structural  Resistance  Term 

Norr  ...  Pressure  j4  f (0) 

Friction  Tenr  ■ f ( sin  0/r.r- 0) 

Equivalent  Fluid  Flow  Term 

Normal  Pressure  ■ f(sin20) 

Shear  Pressure  • f(cos20) 

where  "0"  is  the  local  sJ-'oe  of  the  EP  differential  surface  area. 

The  overall  implication  of  this  effect  is  that  a nose  shape  of  ii  .er- 
mediate  bluntness  may  be  the  best  configuration  for  impacts  Into  hard  targets. 

The  peak  loading  environments  for  this  variation  in  "CRH"  are  presented 
in  Figures  30  through  34.  For  "CRH's"  of  0.5  and  2.25,  twc  peaks  are  given; 
cne  associated  with  chock  (or  compressibility),  and  one  which  muy  be  related 
to  rigid  body  g's.  The  distinction  is  made  because  the  shock  loading  is  a 
short  duration  phenomena  and  would  be  of  interest  to  the  dynamic  structural 
design  of  the  EP,  while  the  rigid  body  loading  would  be  of  nwe  interest  to 
component  inertial  loading.  The  "Shock"  effects  are  only  prevalent  for  the 
low  "CRH"  values  because  of  the  "0"  effect  previously  discussed.  The  penetra- 
tion performance  and  loading  environments  for  the  "blunt"  nosed  EP  are  sum- 
marized in  Figures  35  and  36,  rsspectively.  The  trends  in  penetration  are  as 
expected  with  respect  to  L/D,  which  basically  represents  a weight  change  only 
(i.e.,  normal  impacts  do  not  generate  any  side  loads),  and  performance  is 
somewhat:  below  similar  weight  designs  with  sharper  noses.  Unlike  the  sharper 


nose  design  the  loading  environment  shows  significant  increases  with  velocity. 
This  is  because  of  the  shock  effect  which  is  prevalent  with  blunt  nosed  projec- 
tiles. Accurate  simulation  of  blunt  nose  projectiles  requires  the  consideration 
of  the  effect  of  apparent  mass.  These  phenomena  appear  to  occur  early  in  the 
impact  event  (i.e.,  during  the  initial  phase)  and  consist  of  a wedge  being 
entrapped  on  the  nose  of  blunt  projectiles.  This  wedge  acts  as  an  integral 
part  of  the  EP  and  can  significantly  effect  subsequent  performances.  Refer- 
ring to  Figures  of  the  Series  2-2A,  the  nose  can  be  seen  as  a conical  ex- 
trusion to  the  ’'asic  configuration.  The  3-D  results  for  Runs  2-3,  2-6,  and 
2-18  were  reviewed.  The  peak  surface  tractions  have  been  obtained  from  this 
data  and  are  summarized  in  Table  16.  The  peak  stresses  are  the  same  for  the 
400  and  1000  pound  EP  (as  they  should  be)  with  a fairly  large  increase  in 
stress  going  from  the  1000  to  the  3000  ft/sec  case. 

5.3  SERIES  3 - VARIATIONS  IN  OBLIQUITY  AND  ANGLE  OF  ATTACK 

The  penetration  events  simulated  in  Series  3 were  selected  to  establish 
the  effects  of  obliquity  and  angle  of  attack  on  penetration  performance,  loads 
accelerations,  etc.  The  runs  selected  for  this  series  are  outlined  in 
Table  9.  To  avoid  confusion  and  to  separate  different  effects,  these  runs 
are  grouped  as  shown  in  the  table.  Series  3 consisted  of  Runs  3.1.1  through 
3.64.5  as  defined  in  Table  9. 

The  purpose  of  group  3A  (i.e.,  Runs  3-1  through  3-8)  is  to  determine 
the  effects  of  obliquity  on  the  penetration  parameters.  These  runs  are  sum- 
marized in  Figures  3/  through  39.  The  penetration  depth  (i.e.,  vertical  dis- 
placement) is  also  a designator  of  ricochet  performance.  Although  in  the  runs, 
ricochet  did  not  a'  >ally  occur,  the  shallow  depth  achieved  at  an  obliquity 
of  60  degrees  indicates  that  this  type  of  impact  condition  is  marginal  as  far 
as  ricochet  performance  is  concerned.  The  loading  environments  (rigid  body 
g's)  are  summarized  in  Figure  38.  It  is  apparent  that  the  lateral  loads  in- 
crease smoothly  up  to  60  degrees,  while  only  slight  increases  in  axial  loads 
are  indicated  with  a sudden  drop  off  at  50  degrees.  This  is  because  of  surface 
effects  which  is  shown  in  the  detailed  stress  histories.  Peak  surface  stresses 
remain  pretty  much  the  same  over  this  obliquity  range  and  are  presented  in 
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Figure  39.  The  Integral  of  the  stress  distributions  is  the  parameter  of 
importance. 

Group  3B  (i.e.t  Runs  3.17  through  3.24)  represents  an  angle  of  attack 
effect  investigation  with  some  changes  in  velocity  and  CRH.  This  data  is 
summarized  in  Figures  40  through  44.  The  effect  of  angle  of  attack  on  pene- 
tration performance,  axial  and  lateral  rigid  body  accelerations,  and  surface 
stress  is  shown  in  Figures  40  through  44  for  both  positive  and  negative 
values  of  angle  of  attack,  respectively. 

It  is  apparent  that  penetration  performance  (i.e.  , trajectory  depth)  is 
not  significantly  affected  by  angle  of  attack  as  is  not  the  axial  loading 
environment.  The  lateral  load  however  is  significantly  affected  by  angle  of 
attack  (Figure  42)  and  this  (as  was  verified  by  angle  of  attack  impact  tests) 
has  structural  integrity  implications  relative  to  EF  impact  survival.  A steep 
gradient  occurs  between  a « 0°  and  1°  because  the  angle  of  attack  must  first 
overcome  the  effects  of  obliquity,  which  produce  a lateral  load  opposite  to 
that  of  positive  angles  of  attack. 

Cases  17,  19,  and  20  were  rerun  with  an  increase  in  velocity  to  2500 
ft/sec.  These  runs  are  designated  17A,  19A,  and  20A  and,  referring  to 
Figures  40,  41  and  42  as  expected,  a significant  increase  in  penetration  depth 
was  achieved.  These  same  cases  were  again  run  with  a velocity  of  1500  ft/sec 
and  2500  ft/sec  with  a CRH  of  2.25.  Little  change  in  terradynamic  performance 
occurs  because  of  this  change.  Figures  43  and  44  represent  what  changes  occur 
when  the  angle  of  attack  switches  from  a nose  down  to  nose  up  situation.  For 
this  low  obliquity  situation,  little  change  in  penetration  performance  occurs. 
The  shape  of  the  lateral  accelerator  curve  however  is  somewhat  smoother  than 
for  the  nose  down  angle  of  attack  situation. 

Cases  3.31  through  3.38  involved  a slight  modification  in  the  configura- 
tion and  a change  in  impact  velocity  (i.e.,  1500  to  2500  ft/sec)  as  baseline 
conditions  for  determining  the  effects  of  obliquity.  The  configuration  change 
consisted  of  a small  45  degree  angle  tip  over  a length  of  0.25  inch  on  the 
nose  of  the  EP  instead  of  a pure  ogive  shape  all  the  way  to  the  point.  Runs 
3.31.1  through  3.38.5  are  summarized  in  Figures  45  and  46.  The  overall 
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Figure  43.  Penetration  performance:  negative  attack  angle 
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effect  of  the  tip  by  itself  is  difficult  to  ascertain  because  of  the  signifi- 
cant change  in  velocity,  however  its  overall  effect  is  expected  to  be  umall. 
The  effect  of  this  increase  in  velocity  is  to  delay  the  onset  of  ricochet  by 
approximately  five  degrees. 


Cases  3.39  through  3.47  were  deleted  from  the  original  matrix. 


Cases  3.48  through  3.53  involve  investigation  into  the  performance  of  a 
biconic  EP  configuration.  This  design  is  shown  in  Figure  47.  The  biconlc 
design  represents  a departure  from  the  ogive  nose  type  EP  in  that  resulting 
flow  fields  are  induced  to  separate  at  the  shoulder  of  the  forward  cone  for  a 
period  of  time  during  the  penetration  event.  There  are  two  pertinent  questions 
to  be  asked  and  resolved  relative  to  the  performance  of  this  type  design: 


• For  how  long  a period  during  the  event  is  separation  induced? 


* How  much  of  the  after  flare  remains  in  the  void  generated  by  the 


wake? 


As  a demonstration  of  the  importance  of  these  phenomena,  consider  an  approxi- 
mate loads  comparison  between  the  biconic  design  in  question  and  a design 
which  allows  the  forward  cone  to  continue  out  to  the  full  diameter  of  the  flare. 


Assuming  that  flow  separation  is  induced  at  the  shoulder  of  both  designs  so 
that  the  entire  aft  flare  and/or  cylinder  is  in  the  wake,  the  load  on  the 
cylinder  cone  configuration  is  approximately  proportional  to: 


17 AC  - 


and  on  the  biconic, 


> Prr  Rz 


where 


- resulting  axial  load  on  the  cone  cylinder  design 

■ resulting  axial  load  on  the  biconic 

■ average  normal  pressure 

- cylindrical  radius 


- forward  cone  radius 


2 2 

The  ratio  of  resulting  loads  is  therefore  approximately  R0/RQC  or  in  this 
case,  4.  It  is  apparent  that:  if  separation  did  occur  at  the  shoulder;  and 
if  the  entire  flare  were  within  the  wake;  and  if  the  flare  remained  so  during 
the  entire  penetration  event,  then  the  penetration  performance  would  improve 
by  a factor  of  4.  Actual  performance  is  somewhat  below  this  optimistic  figure. 
At  the  present  time,  there  is  not  sufficient  data  to  define  exactly  the  extent 
of  these  possible  effects.  A limited  number  of  scale  model  test  shots  were 
conducted  into  a hard  target  in  which  an  ogive  nosed  EP  was  compared  to  a 
biconic  configuration  where  the  ratio  R^/R^B  2.4  The  increase  in  penetra- 
tion performance  of  the  biconic  was  only  IS  percent  which  was  u reduction  of 
2,1  less  than  the  area  ratio  factor.  Assuming  that  this  reduction  factor  re- 
mained constant  (a  big  assumption)  then  the  biconic  design  would  have  a per- 
formance increase  on  the  order  of  4/2.1  - 1.9  times  better  than  the  straight 
cone  EP  design.  It  is  Avco's  opinion  that  the  performance  Improvements 
indicated  by  the  above  manner  of  deduction  to  very  optimistic  for  impact 
situations  into  hard  targets.  The  trend,  however,  is  real  and  in  conducting 
Runs  3. 48  through  3.33,  u conservative  assumption  of  the  degree  of  wake  ef- 
fects was  made.  The  summary  curves  for  these  cases  are  presented  in  Figures  48 
through  50.  The  improved  performance  factor  of  this  biconic  design  over  thu 
standard  ogive  EP  computes  to  be  1.08;  i.e.,  Run  1.01  achieved  a penetration 
depth  of  66  inches  while  the  biconic  configuration  went  71  inches. 

The  detailed  trajectory  data  for  Runs  3.60  through  3.64  is  summarized  in 
Figures  51  and  52.  These  runs  were  conducted  to  provide  more  detailed  data 
on  the  effects  of  velocity  and  weight  increases.  As  before,  thu  most  signifi- 
cant trend  is  the  real  improvement  in  penetration  depth  with  velocity  accom- 
panied by  not  too  great  an  increase  in  rigid  body  loads. 

5.4  SERIES  4 - EFFECT  OF  L/D  ON  STABILITY 

Runs  3.17,  18,  20  and  22  were  rerun  in  the  3-D  format  to  provide  surface 
traction  information.  These  data  are  presenced  in  Figures  3.17,1  through 
3.22.10.  The  peak  pressures  are  summarized  in  Table  17. 

The  purpose  of  the  Series  4 parametric  study  is  to  determine  the  effects 
of  EP  L/D  (length  to  diameter  ratio)  on  projectile  stability.  The  target 
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Figure  62.  Axial  and  lateral  acceleration:  velocity  effects 
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madia  in  thin  case  la  a relatively  standard  clay  sand  atx  (glacial  till). 

This  media  has  bees  tested  extensively  at  Avcu  tor  raai*taac*-to-p«o«tratioa 
characteristics.  Consequently , the  results  ot  this  study  can  be  considered 
quite  accurate  and  relatabi*  to  the  real  world. 

The  L/D  range  considered  varied  t ram  4 to  10  wiser*  the  10  vrlue  vm  as- 
sociated with  the  standard  ogive  400  pound  EP.  The  tP  weignt  tor  the  ether 
HP's  was  decreased  proportionally  by  the  L/D  ratio.  Thu  exterior  contigwra- 
tions  of  these  EP's  are  shown  in  Figure  53.  Blunt  nosed  LP's  were  also  con- 
sidered in  this  study  and  are  shown  it  Figure  V4.  The  detail  trajectory  tor 
the  ogive  and  blunt  nose  EP  configurations  is  detailed  in  runs  4.01.1  through 
4.08.5  and  4.09.1  through  4.12.5  respectively  and  sr«  summarised  to  figures 
55  through  59.  It  la  apparent  that  the  most  significant  result  oi  this  study 
is  that  the  low  L/D  EP's  are  not  that  much  less  stable  than  the  higher  L/D 
EP's.  There  are  basically  two  reasons  which  explain  this  result.  They  are: 

• the  center  of  gravity  location; 

• the  penetrated  media  is  till. 

For  EP's,  as  it  is  for  ail  things  that  fly,  the  location  of  the  center 
of  gravity  (eg)  with  respect  to  the  effective  center  of  applied  forces  i»  ot 
prime  importance  to  its  stability.  An  aft  located  eg  will  cause  the  projectile 
to  be  unstable  while  a forward  eg  will  enhance  Its  stability.  During  this 
study,  the  center  of  gravity  was  assumed  to  remain  at  53.5  percent  of  the 
EP's  overall  length.  This  Is  slightly  aft  of  the  mid  point.  When  dealing 
with  actual  EP  designs,  maintaining  this  eg  location  becomes  more  difficult 
as  the  EP  becomes  shorter.  This  is  particularly  true  for  ogive  shaped  EP's 
because  the  ogive  shape  itself  does  not  provide  as  much  volume  up  forward  for 
ballast.  For  this  study  therefore,  placement  of  the  eg  was  somewhat  optimistic. 

The  fact  that  the  penetrated  media  was  till,  aided  in  stabilizing  the 
EP's  because  of  particular  characteristics  relative  to  its  resistance-to- 
penetration.  Till  exhibits  relatively  high  shear  strength  and  friction.  This 
characteristic  results  In  relatively  high  shear  stress  applied  along  the 
length  of  EP  when  the  EP  develops  angle  of  attack.  This  shear  stress  is 
stabilizing  and  contributes  to  a large  degree  to  the  stable  flight  patterns 
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Figure  55.  Penetration  performance:  CRH  = 9.25  and  blunt  nose 
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Figure  58.  Lateral  acceleration:  CRH -9,26 
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of  these  EP's.  This  would  not  be  the  case  in  a media  such  as  air,  water,  or 
dry  sand.  In  such  media,  these  short  EP's  would  probably  be  quite  unstable. 

5.5  SERIES  5 - MEDIA  VARIATIONS 

Series  5 was  reduced  in  scope  from  the  original  matrix  to  include  three 
simulations  into  a reduced  strength  concrete  and  three  simulations  into  till 
(i.e.,  the  same  madia  used  in  Series  4)  and  a prediction  into  welded  tuff. 

The  matrix  of  runs  selected  for  this  study  is  given  in  Table  11.  The  data  are 
summarized  in  Figures  60  and  61.  No  significant  conclusions  were  drawn  from 
this  study  except  that,  in  the  case  of  the  concrete  runs,  this  performance  is 
more  typical  of  what  would  be  expected  in  the  field.  As  discussed  before,  the 
baseline  concrete  resiatance-to-punetration  parameters  used  in  this  study  was 
selected  conservatively  high.  The  value  of  the  trajectory  data  and  the  per- 
formance trends  provided  during  this  study  la  dependent  on  the  credibility  of 
the  tools  used  to  generate  the  information.  As  an  indication  of  accuracy  of 
the  results  a prediction  was  to  be  made  of  an  Impact  event  into  a welded  tuff 
target  at  the  Tonapa  test  site  and  was  to  be  made  prior  to  an  actual  test. 

This  type  of  check  was  done  once  before  into  a soft  target  (i.e.,  the  Watching 
Hill  Test  Site  in  Canada)  with  remarkable  success. * 

The  conditions  of  the  proposed  test  shot  included  the  use  of  the  standard 
ogive  nose  EP,  40U  pounds,  65  inches  long  and  6.5  inches  in  diameter  to  be 
fired  normally  at  1500  ft/sec  into  a welded  tuff  target  media.  The  media 
properties  data  were  provided  by  WES  (Waterways  Experiment  Station)  and  ax*e 
summarized  in  Figure  62.  Avco  fitted  the  bilinear  curve  through  the  data  pro- 
vided and  executed  the  conversion  procedure  described  in  Section  2.0  to  estab- 
lish the  force  law  coefficients  with  which  to  moke  the  required  predictions. 
The  variation  in  target  strength  was  used  in  the  simulations. 

Two  predictions  were  initially  made,  one  at  1500  and  one  at  2000  ft/sec 
with  the  1500  ft/sec  prediction  submitted  to  DNA  prior  to  the  test.  (This  was 
done  in  memo  form  and  is  included  in  Appendix  A,  dated  June  9,  1975.)  The 
impact  velocity  however  was  1640  ft/sec  instead  of  1500  ft/sec  and  the  nose 
of  the  EP  achieved  a depth  of  11.0  feet. 

•Impact  and  Penetration  Study,  Contract  DNA  001-715-C-0181,  Finul  Report  1974, 


Unconfined  compressive  strength  (o^  - Op) 


These  impact  conditions  were  resubmitted  for  a new  estimate  of  predicted 
depth  with  the  result  being  11.3  feet.  The  final  prediction  is  based  on  Runs 
5.1A.1  through  5.3A.5  summarized  in  Figures  63  through  65.  It  should  be 
pointed  out,  referring  to  Figure  65,  that  the  sudden  increase  in  axial  g's  is 
caused  by  the  separation  of  a 157  pound  sabot  used  to  accelerate  the  EP. 

Figures  66  and  67  summarize  the  trend  of  penetration  depths  achieved  within 
the  range  of  1500  to  2000  ft/sec.  It  is  evident  that  the  accuracy  of  the 
simulator  and  the  conversion  techniques  to  establish  the  appropriate  force 
law  coefficients  is  quite  good. 

5.6  SERIES  6 - RICOCHET  PREDICTIONS 

The  purpose  of  Series  6 was  to  establish  the  ricochet  performance  of  two 
operational  bombs,  the  Mk  82  and  Mk  84,  in  a representative  hard  target. 

The  physical  characteristics  of  these  two  warheads  and  their  configurations 
are  shown  in  Figures  68  and  69. 

The  hard  target  is  made  up  of  four  layers  shown  in  Figure  70  along  with 
their  engineering  properties.  Using  the  Cavity  Expansion  Theory  conversion 
procedure,  the  structural  resistance  term  computes  to  be  35,545  psi,  2,536  psi, 
747  psi  and  431  psi,  for  the  four  layers,  respectively.  The  matrix  of  runs 
selected  for  this  study  is  summarized  in  Table  12.  The  impact  velocity  of 
700  ft/sec  was  specified  by  expected  delivery  conditions,  with  the  obliquities 
chosen  to  establif>n  the  ricochet  angle  to  an  accuracy  of  5 degrees. 

The  detailed  trajectory  data  was  generated  in  Runs  6.1A.1  through  6.14c. 3 
and  are  defined  in  Table  12.  The  simulations  which  depicted  ricochet  are 
provided  in  Figures  71  through  76.  Summary  curves  for  the  ricochet  perform- 
ance are  provided  in  Figure  71  for  the  Mk  82  and  Figure  72  for  the  Mk  84.  It 
is  apparent  from  these  curves  that  for  the  0 degree  angle  of  attack  case  the 
Mk  82  experiences  ricochet  at  55  degrees  and  the  Mk  84  at  60  degrees.  The 
difference  is  attributed  to  the  increased  size  of  the  Mk  84  over  the  Mk  82. 
Referring  to  the  summary  curves  (Figures  71  and  72)  a 15  degree  angle  of 
attack  causes  a change  of  approximately  10  degrees  in  the  ricochet  angle.  A 
reasonably  smooth  curve  can  be  estimated  and  drawn  between  the  apparent 
ricochet  points  for  the  angle  of  attack  cases. 
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Referring  to  the  loading  environment  curves,  it  should  be  pointed  out 
that  after  the  onset  of  ricochet  a severe  and  potentially  destructive  loading 
environment  is  applied  to  these  warhead  designs.  This  environment  has  been 
designated  "tail  slap"  and  is  caused  by  excessive  lateral  loads  being  applied 
to  the  tail  section  of  the  warhead  after  the  nose  has  bounced  clear  of  the 
target.  This  type  of  loading  has  been  observed  to  crush  the  back  end  of 
similar  projectiles  and  could  result  in  structural  design  integrity  implica- 
tions relative  to  the  Mk  82  and  84.  In  the  case  of  aft  mounted  fuzing  systems 
this  could  cause  problems  relative  to  fuze  functioning. 
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APPENDIX  A 


9 June  1975 
C500-75-DH-53 


Major  Todd  Stong 
Defense  Nuclear  Agency 
Washington,  D.C. 

Dear  Todd, 

Enclosed  you  will  find  Avco's  terradynamic  performance  predictions 
of  the  1500  ft/sec  Tonapa  DEP  Test.  The  structural  reals tance-to- 
penetratlon  parameter  " »/  " was  based  on  a conversion  process  using 
the  Cavity  Expansion  Theory  with  the  basic  data  being  supplied  by 
WES.  The  inertial  resistance-to-penetration  was  based  on  our  own 
estimates  of  the  characteristics  of  Welded  Tuff. 

Referring  to  Figure  A-3,  the  estimated  penetration  depth  is  a little  over 
9 feet. 

Accelerations  are  obtained  from  Figure  A-4.  Initial  peak  acceleration  is 
on  the  order  of  3700  g while  the  launch  sabot  is  still  attached.  The 
slope  to  4600  g is  caused  by  the  assumption  that  the  sabot  does  not 
detach  instantaneously. 

Relative  to  the  ongoing  "Impact  and  Penetration  Technology  Program," 

Teats  3 and  4 were  conducted  this  week  (i.o.,  June  27th  and  28th)  and 
they  both  went  well. 

Teat  3,  which  consisted  of  an  impact  velocity  of  1500  ft/sec  and  an 
angle  of  attack  («  ) of  2 degrees  produced  stresses  (based  on  strain 
gage  readings)  as  expected  in  between  those  of  u » 0°  and  5°.  Test  4, 
which  consisted  of  an  impact  velocity  of  1500  ft/sec  and  a «■  10°  produced 
peak  structural  stresses  of  approximately  300,000  psi,  again  a nearly 
linear  extrapolation  of  the  other  tests. 
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Figure  A-2.  Tonapa  rest  prediction  trajectory  data 
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